Averaged collision and reaction rates in a two-species gas of ultracold fermions 



(N 

o 

(N 

> 
O 

oo 



C/2 ■ 

i 

^— > ■ 

a- 



o 
o 



X 



Alexander Pikovski 

Institut fur Theoretische Physik, Leibniz Universitat Hannover, Appelstr. 2, 30167 Hannover, Germany 

Reactive or elastic two-body collisions in an ultracold gas are affected by quantum statistics. 
In this paper, we study ensemble-averaged collision rates for a two-species fermionic gas. The two 
species may have different masses, densities and temperatures. We investigate how averaged collision 
rates are affected by the presence of Fermi spheres in the initial states (Pauli blocking of final states 
is not considered). It is shown that, independently on the details of the collision, Fermi- averaged 
collision rates deviate from Boltzmann-averaged ones, particularly for a gas with strong imbalance 
of masses or densities. 



I. INTRODUCTION 

Experiments with ultracold atoms and molecules probe 
collisional physics at very low temperatures. Gases of 
quantum degenerate fermionic atoms are being stud- 
ied |l| , and experiments are progressing towards the cre- 
ation of degenerate gases of fermionic molecules Q. At 
very low temperatures, quantum statistics plays an im- 
portant role in interparticle collisions. One way to ob- 
serve an effect which depends on the symmetry of the 
particles is to look at low-energy collisions which obey 
quantum threshold laws; this was measured in elastic 
collisions in fermionic atomic gases and in reactive 
collisions of ultracold molecules Q. Another effect of 
quantum statistics is Pauli blocking of collisions, as was 
observed in Ref . [l2j : collisions are suppressed if the final 
state cannot enter a filled Fermi sphere. In contrast to 
the quantum threshold behavior, the observation of this 
effect requires the fermions to be quantum degenerate. 

In the present work we investigate how averaged col- 
lisional rates depend on the temperature, density, and 
masses in a two-species Fermi gas. The usual experimen- 
tal observable in a cloud of ultracold gas is a collision rate 
which is averaged over the distribution of velocities in the 
gas. Usually this ensemble-averaging of collision rates is 
done using a Boltzmann distribution, however as experi- 
ments move deeper into the quantum degenerate regime, 
quantum statistics starts playing a role. Particular atten- 
tion here will be focused on the imbalanced Fermi gas, 
i.e. a two-component gas where the two species have dif- 
ferent densities or masses. The effect of mass or density 
imbalance on the averaged collision rates at low temper- 
atures is most pronounced in this setting, the rates can 
differ significantly from a Boltzmann-averaged rate. 

Averaged collision rates for a two-component gas of 
fermions are studied, for general two-body collisions with 
an arbitrary cross-section. The particles are assumed 
to be in equilibrium; effects of Pauli-blocking are not 
considered. The two types of fermions may have dif- 
ferent masses, different densities, and different tempera- 
tures. First, it is shown that an averaged collision rate 
is expressed as an average over the distribution of rela- 
tive energies (SecQl]). This distribution is calculated for 
fermions analytically at T = 0, and it is shown that this 
function has a simple geometric interpretation fSec. Mil) . 



At finite temperatures, the relative energy distribution 
is calculated numerically and a simple approximation for 
all temperatures is proposed (Sec. IHIj) . In Sec. [V] we dis- 
cuss the threshold laws for ensemble-averaged collision 
rates. Some details of the calculations are presented in 
the Appendices. In Sec. IIV1 we discuss how the results 
can be transferred to a situation where the particles are 
confined to a two-dimensional geometry. 



II. AVERAGED COLLISION RATE 

Experiments with clouds of ultracold gases measure 
collision rates, e.g. using photoassociation spectra [l3| or 
trap loss rates due to chemical reactions Q . Here one ac- 
tually measures rates which are averaged over the distri- 
bution of velocities in the gas. This ensemble-averaging 
is the topic of this section. 

Consider a collision between two particles, A and B, 
which results in some final state /: 



A + B -> /. 



(1) 



This may be an elastic collision, an inelastic collision, 
or a reaction. The particles A, B have masses rriA, tub 
and velocities va, vb in the laboratory frame. The total 
scattering cross-section for the process (JXJ) in the center- 
of-mass frame will be denoted cr(p), it depends on the 
relative momentum p of the colliding particles. 

In a gas cloud, where many collisions take place, the 
number density n a of particles in state a = /, A, B sat- 
isfies the rate equation 



dn f i i w 

—— = (va[p)) n A n B 

at 



(2) 



where (va) is the two-particle collision rate averaged over 
the distribution of velocities in the system. The rate 
equation ^ refers only to the process (p}. It is assumed 
that, in the case of fermions, the final state / is not oc- 
cupied, thus there is no effect of Pauli blocking. If the 
particles A and B are indistinguishable, a small modifi- 
cation to Eq. © may be necessary, see @. 

The relevant quantity for many experiments is the 
ensemble-averaged collision rate K = (va), it has dimen- 
sions (length) 3 /time. If velocities of particles of type A, 
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B are distributed according to the distributions f A , f B , 
the averaged collision rate K is @ 



K = (v<r(p)} 



va(p)f A (v A )f B (v B )d 3 v A d 3 v B (3) 



We will consider the case where the velocity distributions 
and the cross-section depend only on the magnitude of 
the velocities, i.e. they depend only on the energy of the 
particles. Changing to the center-of-mass frame with the 
formulas of Appendix [A] brings the integral to the follow- 
ing form: 



K=J- f E Icl a(E rcl )F{E rcl ) dE Icl 
V m J 



(4) 



with 



1 pit POO 

F{E, cl ) = - / MsmO f A {E A )f B {E B )E^dE cm . (5) 

2 Jo Jo 

Here f AjB are the energy distribution functions. The 
normalization for the velocity distribution functions 
is Jf(v)d 3 v = 1, the energy distribution functions 
are normalized to / °° f(E)^/EdE = 1, and f A _ B — 
2-K{2jm A ^ B ) 3 l 2 ] A ^ B . The expressions for E A and E B in 
terms of E cm , E Te \ and 8, which are needed to evaluate 
the integral (|5|), are given in Eq. (|A3I) . 

The ensemble-averaged collision rate is given by the 
cross-section a averaged with the function F, as shown 
by Eqs. g]), ©. The function F(E m i), which we call 
distribution of relative energies, contains all the depen- 
dence on the densities and temperatures of the particles 
involved in the collision. The cross-section cr(p) contains 
all the information about the two-body scattering. 

The expression for the distribution function of rela- 
tive energies, Eq. ([5]), can be brought to a different form 
which will be of use later. Transforming the integral in 
Eq. ([5]) to the coordinate system (x,y, z), as described in 
Appendix [AJ we have the representation 



F{E m} ) = 



f A (E A )f B (E B )dxdydz. (6) 



The integration runs over the whole space, E A and E B 
are expressed through x,y,z using Eqs. (IA3|) . (IA5|) . and 
the mass factor r/ is defined in Eq. (|A"6j) . 



III. DISTRIBUTION OF RELATIVE ENERGIES 
FOR FERMIONS 

Now we discuss the distribution function of relative 
energies F(E TC \), which determines the averaged colli- 
sion rate [cf. Eq. flU)] for a two-species gas of fermions 
in equilibrium. We will consider a Fermi gas where the 
two species may have different densities or masses. Such 
imbalanced systems have attracted recent experimental 
11911 and theoretical !I8J attention. 




FIG. 1: Intersection of two spheres with radii r\ and T2, cen- 
ters separated by distance d. 



A. Zero temperature 

First, the case of zero temperature is considered. The 
distribution functions for the particles A and B are nor- 
malized T = Fermi functions 



HE) 



2^/2 



(7) 



where 9{x) is the step function and /i is the chemical 
potential. The particles of type A, B may have differ- 
ent chemical potentials fj, A , fi B . It will be convenient to 
introduce the "reduced" chemical potentials 



v A 



2m B 



m A + m B 



-HA, 



2m A 



m A + m B 



Hb, 



(8) 



and for equal masses we have v A ^ B = (i A , B . When the 
Fermi functions, Eq. ([7]), are inserted in the represen- 
tation ([B]) for the distribution of relative energies, one 
encounters the volume integral 



9{ha - E A )9(]i B - E B )d 3 x. 



(9) 



Once E A and E B are expressed in terms of x,y,z, it seen 
that this integral has the following geometric meaning: 
it is the volume of intersection of two spheres. Let us 
briefly discuss this volume. 

The volume of intersection of two spheres, which have 
radii r\ and and whose centers are separated by a 
distance d (see Fig. [T]), is given by the expression 



V L (d; r u r 2 ) = —(r 1 +r 2 -d) 2 {d 2 +2d(r 1 +r 2 )-3(r 1 ~r 2 ) 2 }. 

This formula, as it stands, is only valid as long as the 
volume of intersection has the shape of an asymmetric 
lens, as indicated in Fig. [TJ However, if the distance 
between the spheres d is large, there is no intersection, 
and if d is small, one sphere is located completely inside 
of the other. We can write for the volume of intersection 



V{d-r u r 2 ) 







,3 4 



mm(g7rrf 
V L {d\ri,r 2 ) 



if d 2 > (n +r 2 ) 2 
7rr|) if d 2 < (ri - r 2 ) 2 
otherwise 

(11) 
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FIG. 2: F(E rel ) for a) v A = 0.25, v B = 1.0 and b) v A = vb = 



1. The dashed lines are at V ^+" B 
y/uAfB, the dotted line is at /i» = — 



\Jvavb and 



"A+ V E 



+ 



A + V B 



In this form the expression is valid for all values of d, r% , 
and T2. 

Returning to the distribution function of relative en- 
ergies, it is seen by using ([7]) in Eq. (j6|) that the result 



rei; \VA, \JVB 



(12) 

This is the volume of intersection of two spheres with 
radii \pUA an d \Jvb , the centers of the spheres being sep- 
arated by y/2E ie i, up to normalization. It is interesting 
to note that the masses of the particles do not appear 
explicitly, but enter only through the reduced chemical 
potentials i>a, vb- 

Let us discuss the general form of the function Fo(E Te \). 
A plot is shown in Fig. [2^ for the strongly imbalanced 
case v\ <C vi- This case occurs either if the particles have 
different masses or if the particle densities are different, 
cf. Eq. ©. The function has the form of a smoothed 
step. From the cases in Eq. (1TT1) . it is seen that Fo(E rc \) 
is a decreasing function in the range (shown by dashed 
lines in Fig. [2]) 



va + v B 



- Vvavb < E rcl < 



+ Vv A v B - (13) 



For smaller values of E rc \ the function is constant and 
for larger values of E le \ it is identically zero. The step is 
centered at /i* = (y A + ^s)/2, but it is not symmetric 
around this point. 

In the balanced case [i = v A — vb, the expression for 
-Po^rci) simplifies considerably. This case is shown in 
Fig. [5Jd. The distribution function has no flat region, it 




T=0 

T=0.2 

T=0.5 

T=0.9 

T=1.5 

T=1.5, Boltzmann 



FIG. 3: F(E le i) at different temperatures, for ua = 0.6, vb 
1.5 at T = 0. 



is now a strictly decreasing function between E le i = 
and E re \ = 2/i. 



B. Finite temperatures. 

For finite temperatures, the distribution functions of 
the particles are normalized Fermi-Dirac distributions 



fT(E)=Af^T 



1 



1 + e (E- M )/T' 



(14) 



the normalization constant AT is given in Appendix 151 
We allow the particles A, B to be at different temper- 
atures Ta-, Tb- The corresponding distribution of rel- 
ative energies Ft(E ic i) becomes, geometrically, the vol- 
ume of overlap of two spherical distributions whose radii 
are broadened by temperature. One can introduce the 
reduced temperatures 



2m 



t A 



B 



m A + TUB 



-Ta, 



tb 



2niA 



m A + m B 



-TV 



(15) 



As shown in Appendix [B] the distribution of relative en- 
ergies can be written as 



F T (E le 



Af VA , ta Nv B , t b 



36 TATB 

{e A e B ) 3/2 F (E rcl ;€ A ,€B) 



, -, ^de A des- (16) 

o cosh 2 (^)cosh 2 (^) 



Note that the masses of the particles enter only through 
the reduced chemical potentials v a and temperatures r a . 

In the limit of high temperatures, both particles are 
distributed according to the Boltzmann distribution 



foo(E) = 



,-E/T 



(17) 



In this case the distribution of relative energies may be 
evaluated directly from Eq. ([S]). Then we obtain, as 
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known Q, that the distribution of relative energies is 
a Boltzmann distribution 



Foo(E Te i) — 



3/2 



with the effective temperature 



T» = 



ta + t b _ m A T B + m B T A 
2 mA + rriB 



(18) 



(19) 



The distribution function of relative energies may be 
calculated numerically from Eq. ([5]) or from Eq. (fl6|); for 
low temperatures the form in Eq. (|16|) is more conve- 
nient. Figure [3] shows the result of such a calculation. 
The chemical potentials va 1 v B were fixed at T = 0, 
and at T > the temperature-dependent chemical po- 
tentials were found by solving Eq. (|B5[) . The figure also 
compares the distribution function with the Boltzmann 
[-FooC^rci) from Eq. JTSJ] limit. It is seen that the dis- 
tribution function approaches the Boltzmann form when 



C. Approximation for all T 

Since the distribution of relative energies has a com- 
plicated form at finite temperatures, an approximation 
might be helpful. The distribution function of relative 
energies may be approximated, at all temperatures, by a 
Fermi distribution 



F 



T,app 



r(-E'rel) = -^A*.. 



1 



1 + e (s„i-M»)/r» 



(20) 



with the effective temperature T* from Eq. 
effective chemical potential 



T9l and the 



va + v B 



tuaPb + tub Ha 
m A + m B 



(21) 



This gives a reasonable approximation for the shape 
of Ft(E ic \) at low temperatures. For high tempera- 
tures, it approaches the exact distribution F 00 (E Ie i) from 
Eq. (IB). 



IV. TWO DIMENSIONS 

Gases of ultracold atoms or molecules can be con- 
fined to thin layers. For example, recent experi- 
ments at JILA |8j investigated a gas of KRb molecules 
in a pancake-shaped geometry and measured the rate 
of chemical reactions; related theoretical work consid- 
ered ensemble-averaged rates (Boltzmann distribution in 
Ref. [HI and Fermi-Dirac distribution in Ref. [13]) for a 
two-dimensional geometry. When the gas is confined to 
a thin layer, one can describe collisions as if the parti- 
cles were moving in a two-dimensional plane. The col- 
lisional cross-section, however, depends on the details 



of the interparticle interaction and the confinement. In 
this Section we will consider how the preceding results 
(Sec. HI and Sec. Ifflj) can be adapted to the case of two- 
dimensional scattering. 

For two dimensions, what was said at the beginning 
of Sec. [U remains valid, but now in Eq. ([2]) one should 
use the two-dimensional number densities (dimensions 
1 /length 2 ), and the two-dimensional scattering cross- 
section <t' (dimensions 1/length). We will add a prime to 
two-dimensional quantities to avoid confusion. 

The averaged two-dimensional collision rate K' , having 
dimensions (length) 2 /time, is 

K' = {va'(p))= ( ' [va'(p)f' A (v A )f' B (v B )d 2 v A d 2 v B . (22) 



For velocity distributions which depend only on the en- 
ergy of the particles, changing to the center-of-mass (see 
Appendix [3} brings the integral to the following form: 



/ Q /-OO 



with 



F'(E iel ) = j — f A (E A )f' B (E B )dE c 



(23) 



(24) 



Here f' A B are the energy distribution functions, the nor- 
malizations are Jf'(v)d 2 v = 1 and J °° f'(E)dE = 1, so 

f' AB — (2Tr/nT>A,B)fA b ■ T ne expressions for Ea and E B 
in terms of E cm , E Te \ and 9, which are needed to evaluate 
the integral (|2~4"|) . are given in Eq. (|A3I) . 

The double integral in Eq. ([M]) can be transformed 
to the coordinate system (x 1 , y'), see Appendix [Al which 
results in 



F'(^ re i) = ^JJ f A (EA)f B (E B )dx'dy'. 



(25) 



The integration runs over the whole space, Ea and E B 
are expressed through x',y' using Eqs. (|A3|) and (|A8jl . 
and rj is given in Eq. (|A6|) . 

So far the discussion has been quite general, now we 
consider the case where the particles A, B are fermions. 
At zero temperature the distribution functions for the 
particles A, B are 



f' {E) = -6(n-E), 



(26) 



with possibly different chemical potentials fiAi Hb- The 
integral that results when (121)1) is inserted in Eq. (|2"5j) has 
the geometric interpretation of the area of overlap of two 
circles. For finite temperatures, the boundary of these 
circles becomes broadened by temperature, as before. 

The area of intersection of two circles, with radii r\ 
and T2 separated by distance d (see Fig.rjJ), is 



.4, 



-Ti arccos 



2dri 



-r 9 arccos 



2dr 2 



\^(rl+rl+d?f -2(rt + rl+d*). (27) 
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To include the cases where the circles are completely dis- 
joint and where one circle is inside the other, we write 

ro if d 2 > (n + r 2 ) 2 

A(d;ri,r 2 )=< min(7rr 2 , Tirf) if d 2 < (n - r 2 ) 2 
I, Ax,(d;ri,r2) otherwise 

It follows that the distribution function of relative en- 
ergies at T = is 

2 



rcl J 



(29) 



where the v a are given in Eq. ([5]). This function has 
a similar shape as the corresponding function in three 
dimensions. 

In the high-temperature case, the energies of the par- 
ticles A, B follow the Boltzmann distribution 



fL(E) 



1 

T ' 



-E/T 



(30) 



with possibly different temperatures Ta, Tb- Then the 
distribution of relative energies, calculated from Eq. (|2"4"1) . 
becomes a Boltzmann distribution in the relative energy 



-E IBl /T, 



(31) 



with the effective temperature T* as in Eq. (fT9|) . This 
result has the same form as in three dimensions. 



V. THRESHOLD LAWS 

The dependence of the averaged collision rate on the 
parameters of the two-species gas (mass, temperature, 
chemical potential) can be found using the results of 
the preceding sections, if the scattering cross-section 
is known. Some conclusions can be drawn, however, 
without specific knowledge of the cross-section in the 
low-energy collisional regime where quantum threshold 
laws a pply . According to the Wigner threshold laws 
(see [illlCll] for a review) , the cross-section at low energies 
scales with the relative energy as 



a(E) oc E p , 



(32) 



where the exponent p depends on the type of collision. 
For elastic collisions with short-range forces in the £-th 
partial wave channel, we have p = 21, If molecules un- 
dergo chemical reactions upon collisions, the quantum 
threshold laws hold only if the chemical reaction barrier 
is absent or much lower than the centrifugal barrier of 
scattering 0, For these exothermic barrier less re- 

actions in the t-th. partial wave, p = £ — 1/2. In gases 
of ultracold atoms or molecules, one expects collisions 
to take place in the lowest partial wave channel that is 
allowed by exchange symmetry. 

The averaged rate constant at zero temperatures for 
collisions obeying Eq. (|3"2")l can be obtained in closed form 



from Eqs. ([4]) and (fl~2|) . Here we only note the following 
property of this rate Kq. If the chemical potentials are 
scaled as va — > Az/a, vb — !> ^vb, then Kq \ p+1 / 2 Ko- 
In particular, for p = —1/2 this means that Kq is con- 
stant. 

For finite temperatures, the exact dependence of the 
averaged rate constant on the reduced temperatures ta, 
tb and chemical potentials va, vb can be determined 
from the representation (fT6]) . Inserting Eqs. (|16l) and (f32|) 
into Eq. ^ and performing the integration over E TO \ first, 
one obtains expressions for K. For p = —1/2, one finds 
again that K is constant, and for other values of p more 
complicated expressions result which involve combina- 
tions of Fermi-Dirac integrals (see Appendix [B] for their 
properties). In the strongly quantum degenerate regime 
(/i/T 3> 1), one can obtain a low-temperature expansion 
applying Eq. flB4|). 

A qualitative picture of the dependence of the averaged 
rate constant on the chemical potentials and temperature 
can be obtained by approximating F{E rc \) with a Fermi 
function with chemical potential [i* and temperature T* , 
cf. Eq. (|2U)l . Then, the averaged collision rate for cross- 
sections obeying Eq. (|3"2j) has a simple form, it scales for 
all temperatures as 



-^appr (X T* 



p+1 



(33) 



where the T are the Fermi-Dirac integrals. 

For high temperatures, the distribution of relative en- 
ergies approaches a Boltzmann distribution. The thresh- 
old law (|3"!2")) then leads to the following well-known scal- 
ing of the collision rate 



KotT p+1/2 , 
which also follows from Eq. (f33|) . 



(34) 



VI. CONCLUSIONS 

Two-body collision or reaction rates in an ultracold 
gas are affected by quantum statistics. The partial wave 
channel of the collision depends on the symmetry of the 
particles involved; the density of final states may be mod- 
ified (e.g. Pauli blocking for fermions); and, finally, the 
available initial states depend on the ensemble. In this 
paper, we have investigated initial-state effects for a gas 
of two-species fermions with possibly different masses, 
particle densities, and temperatures. Independently of 
the details of the collision, the ensemble-averaged colli- 
sion rate is the cross-section averaged with the distribu- 
tion of relative energies. This distribution function was 
calculated and its form was discussed in detail. An ap- 
proximation was proposed, which becomes exact in the 
high-temperature limit. It was also shown how to trans- 
fer these results to particles confined to two dimensions. 
As an application of the results, we considered the scal- 
ing of the averaged collision rates for the case where the 
cross-section follows the Wigner threshold law. 
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The results will be of use for the analysis of experi- 
ments which measure collision or reaction rates in gases 
of ultracold atoms or molecules. In particular, they will 
serve to determine how the measured averaged rates in 
quantum degenerate two-species gases depend on the 
temperature, the masses and densities of the particles. 

I would like to acknowledge helpful discussions with 
L. Santos, M. Klawunn and C. Salomon. 



Appendix A: Transformation of coordinates 

The coordinate transformation between the laboratory 
frame and the center-of-mass frame is given by the fol- 
lowing expressions. Two particles A, B have masses tua, 
m B and velocities v^, vb in the laboratory frame. The 
velocities in the center-of-mass frame are 

v = va-vb, {jiiA + m B )V = vtlava + m B v B , (Al) 

The corresponding momenta are p = m(yA — v b) and 
P = (rriA + ms)V, where m = m^rriB / (tua + ms). The 
related energies are 



TP 1 2 



w 1 2 

E B = 2 m BV B > 



(A2) 



E, 



cm - A m A + m B )V , E rel = -mv 



The angle between v and V is denoted by 9. The en- 
ergies in the laboratory frame and center-of-mass frame 
are related by 



m vn 

Ea = ficm H E xe \ 

m B rn,A 



4mE rr i E. 



rcl-^cm 



rriA + m B 



cos 9 



mm AmE rc \E cm 

E B = E cm + E Icl - W ■ cos 9 

mA m B V mA + m B 



and also 



E r 



E re \ = Ea + El 



(A3) 



(A4) 



The Jacobian of the transformation from (va,Vb) to 
(v, V) is unity. 

To discuss the integral over the center-of-mass energy 
E cm and the angle 9 in Eq. ([5]) (here E le \ is a regarded as 
a constant parameter) , we view 9 as the azimuthal angle 
in spherical coordinates, and introduce the coordinates 
(x,y,z) as 



x 



= vVe~ c 



z = rj\J E cm cos 9, 
with 



The Jacobian is 



3in6>, y = <q^/E c 



V / 2mAm B 
mA + m B 



sin m sin f 



(A5) 



(A6) 



d{x,y,z) 1 3 . r— 
m ^ ¥) = - 2 V *aOVEZ. (AT) 



The coordinates (x, y, z) are proportional to the Carte- 
sian components of the center-of-mass velocity V or the 
center-of-mass momentum P. 

For the two-dimensional integral over the center-of- 
mass energy, one can use the coordinates (x',y') given 
by 



x' = r}yjE cm cos 6, y' = rj \J E cnl sin 9 (A8) 
with the Jacobian 



d(x',y') _\_^ 



d(E cm ,0) 2 



(A9) 



Appendix B: Distribution function at finite 
temperature 

The distribution function of the particles in a Fermi 
gas is the normalized Fermi-Dirac distribution 



f T (E)=Af ti ; 



1 



I + e (E-/x)/T ' 

with the normalization constant 

2 



such that J™f(E)VEdE = 1. Here 

1 f°° P 

Fj{x) 



T{l+j)J 1 + e 



-dt 



(Bl) 



(B2) 



(B3) 



is the Fermi-Dirac integral of order j [16l | - The expan- 
sion of this function for the strongly degenerate, low- 
temperature regime (x 3> 1) has the form 



Tj(x) 



rJ+l 



r(j 



[l + C: 



(B4) 



for half- integer j, see Ref. [16[ for more details. For the 
non-degenerate, high-temperature regime (x <C —1) we 
have Tj(x) = e x . 

The chemical potential of a Fermi gas, for a fixed par- 
ticle density, depends on the temperature. It is related 
to the Fermi energy Ep = /i(T = Q) by 



Oi/T) 



4 /E F \3/ 2 



30FV T 



(B5) 



In order to derive a representation for the integral ([5]) 
when j a and f B are Fermi-Dirac distributions, write for 
the Fermi function np(£) = (1 + e^/ T ) _1 : 



M0 = -/' 



— {u)du = - 9(u-£) — (u)du 



Here 



driF 



-IT cosh 2 f — 
2T 



(B6) 
(B7) 
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Inserting Eq. (IB6|) into Eq. ([5]) and changing the order of of relative energies which was discussed in Sec Mil The 
integration, we encounter the T = distribution function result is given in Eq. (|16p . 
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